sra (version 0.1.4.1)

6. sraData: Generates an object of class "sradata", necessary to run the models provided by the sra package.

Description

The data necessary to analyse selection response time series are, for each generation, (i) the mean phenotype of the population, (ii) the phenotypic variance, (iii) the mean of the breeders, and (iv) the population size. These data have to be grouped into an object of class "sradata", which can be provided to the sra analysis functions such as sraAutoreg or sraCstvar.

Usage

sraData(phen.mean, phen.var, phen.sel, var.sel=NULL, N=NULL, 
  gen=NULL, rep=NULL)

Value

An object of class sradata.

Arguments

phen.mean

The vector of phenotypic means.

phen.var

The vector of phenotypic variances.

phen.sel

The vector of the mean phenotype of breeders. Can be NA for the last generations.

var.sel

The vector of the phenotypic variances if breeders. If not provided (NULL), the vector is estimated assuming truncation selection and normal distribution of phenotypes.

N

The vector of population size for each generation (before selection).

gen

The generation numbers. Useful when several times series are provided.

rep

The repetition identification. Useful when several time series are provided.

Details

If not provided, the default value for N is 100. Incorrect values for N will affect the likelihood value and the maximum-likelihood estimates.

gen will be assumed to vary from 1 to the maximum number of generations by default.

If more than one time series are provided, it is dangerous not to specify gen or rep. rep can be any unique identifier. For instance, fro two times series of 3 generations each, gen can be c(1,2,3,1,2,3) and rep can be c("up","up","up","down","down","down").

See Also

sraAutoreg, sraCstvar, and other mechanistic models, sraAutoregTimeseries.

Examples

########### Generating a dummy dataset ################

m <- c(12,11,12,14,18,17,19,22,20,19)
v <- c(53,47,97,155,150,102,65,144,179,126)
s <- c(15,14,14,17,21,20,22,25,24,NA)
n <- c(100,80,120,60,100,90,110,80,60,100)

########## Making a sra data set #######################
data <- sraData(phen.mean=m, phen.var=v, phen.sel=s, N=n)